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Functional sensitivity of Burgers and related equations to initial conditions 
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In this paper, we apply sensitivity methods to nonlinear PDEs like Burgers and KPZ equations. These equa- 
tions are known to have analytical solutions which make easier the analysis of the sensitivity of their solutions 
to initial conditions. The main result stands in the fact that the most the solution is sensitive to the initial condi- 
tion, the most it is decorrelated in space, i.e. the values of the initial condition participate to the solution at all 
distances of the wave front. This finally reveals a particular aspect of the Burgers turbulence. 



I. INTRODUCTION 

Burgers and related equations have been introduced in 
many fields of sciences such as non-equilibrium statistical 
physics. For instance, we can cite, in cosmology, the model 
known as the adhesion model 1 1], there is also the modeling 
of traffic jam |2|. the description of directed polymers in ran- 
dom media or the dynamics of growing interface |5]. 
Note also that the Burgers equation may be a draft modeling 
to fluid dynamics. 

When modeling a system by means of the Burgers or the 
KPZ equation knowledge of the initial condition is required. 
However a great insight of the problem is necessary in or- 
der to fit experimental data as soon as high quality results are 
available. It is then interesting to know the interdependence 
of the solution of the modeling equation to initial conditions. 
This is the goal of sensitivity analysis to provide the system re- 
sponse to variations of input |6, 7]. This mathematical method 
has been applied in numerous domains of sciences, see for in- 
stance 

For an initial value problem, such as the following evolution 
equation: 



dtu(x, t) = A u(x, t) 
u(x, 0) = <p(x), x El 



(1) 



where A is some nonlinear operator, the variation 6<p(x) of 
the initial condition should imply variations Su(x, t) in the 
solution u(x,t) of the equation. These variations obey the 
functional relation: 



5u(x, t) 



Su(x, i) 
6<p(y) 



Stp(y)dy. 



(2) 



Then the functional derivative 5u(x,t)/§cp(y) (sometimes 
called the density) gives a quantitative measure of the re- 
sponse of the actual solution u(x,t) to any variation of the 
input. 

As said extensively in the literature, nonlinear and chaotic 
systems are mainly characterized by sensitivity to initial con- 
ditions. The purpose of the present paper is to quantify this 
sensitivity in the cases of Burgers and related equations. 



II. THE HEAT EQUATION 

Before we treat the case of nonlinear equation, we shall 
examine in this section the case of the heat equation: 



d t Z = vd xx Z, 



(3) 



with the initial condition: Z(x,0) = i/j(x), and where v is 
some diffusion parameter. We have the well known solution 
to the heat equation as the convolution integral: 



Z(x,t) = 



1 



\J\tuA 



ip(y) exp 



(x 



dy. (4) 



Now we can calculate the sensitivity coefficient, i.e. the func- 
tional derivative of the solution Z(x,t) with respect to the ini- 
tial condition ip(x) directly from this solution. It reads: 
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6Z(x, t) 



H{y) 
5i){x') 



exp 



(x - y) 2 



Avt 



Ay. 



As we have 
then given by: 

5Z(x,t) 



(5) 

5(x' — y), the sensitivity coefficient is 



Sijj(x') y/Anvt 



exp 



{x-x'f 



(6) 



We first remark that the density is well correlated along the 
line x = x' in the (x, x') plane and it does not involve the 
initial condition ip(x) itself. In fact, since the heat equation 
is linear, this density is also a solution to the heat equation 
with the initial condition S(x — x') {i.e. the fundamental so- 
lution). Moreover, as t goes to infinity the sensitivity coeffi- 
cient spreads and goes to zero. Consequently, the solution to 
the heat equation is asymptotically insensitive to initial con- 
ditions. This proves, as well, the known issue concerning the 
very high difficulty to find the initial condition from the mea- 
surement of Z(x, t) at any time. 



III. THE UNFORCED BURGERS EQUATION 

Let us recall, the chain of transformations leading to the so- 
lution of the unforced Burgers equation. The Burgers equation 
is given by: 



d t u + ud x u = vd xx u, 



(7) 



2 



with the initial condition u(x,0) — ip(x). Where v stands 
usually for the viscosity coefficient. The change of function 
defined by u(x.t) = —d x h(x,t) leads to the KPZ equation 



(8) 



with the corresponding initial condition h(x, 0) = rj(x). Then 
the Hopf-Cole transformation: h = 2v In Z, yields to the heat 
equation (Eq. ®) lH El • 

From the general solution of the heat equation (Eq. 0), we 
find the solution to Eq.© as: 



h(x,t) = 2i/ln 
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2t 
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dy } . (9) 



Taking the functional derivative of this equation, we obtain 
the sensitivity coefficient of h(x, t) to the initial condition: 



5h(x, t) 
5r/(x') 



exp 



It 



dy 



(10) 



On the contrary to the heat equation, the unforced KPZ 
equation solution is very sensitive to the initial condition, 
moreover this is also the case in the long time limit: 



lim 

t — >oo 



exp 



8h{x, t) 

~^M = JL^exp^)] dy 



(11) 



where we assume the convergence of the integral. Surpris- 
ingly it does not depend on x. This means that, in the long 
time limit, the initial condition r/(x') at the distance x' influ- 
ences h(x, t) for all the values of x with the same weight. It 
is interesting also to look at the inviscid limit of this result. 
Assuming there is only one stationary point a, the sensitivity 
coefficient may written as: 



6h(x 7 1) 



5rj(x') 



h"(a) \ 
iiru 



■ exp 



2// 



(r)(x') - r](a)) 



(12) 

Note these two limits do not act always uniformly. As a cor- 
responds to the maximum of rj(x), we see that the sensitivity 
coefficient is peaked on this value. 

A simple example of illustration of these results may be 
given with the choice of the initial condition for the velocity 
in the Burgers equation as u(x, 0) = sech 2 x, so that the ini- 
tial condition for h is h(x,0) = — tanh x. On Fig. (Q- a 
and b), we present the values of h(x,t) and of the solution 
to the Burgers equation u(x,t) for the time t = 8 (given in 
arbitrary units) as functions of x. The dotted lines correspond 
to initial values of these functions. The viscosity is fixed to 
the rather small value: v = 0.05. The sensitivity coefficient 
8h(x, t)/Sr](x') as a function of x is given on Fig. c) for 
a fixed value of x' = 4. We observe a maximum value of the 
sensitivity at the bottom of the wave front. This is also the 
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FIG. 1: Plots of the solution to Eq. j5J: (a) and of the solution to the 
Burgers equation (b) for the time t = 8 in arbitrary units. The doted 
lines correspond to initial conditions. The sensitivity coefficients are 
plotted as function of x on (c) and as function of x' on (d). 



case for the sensitivity coefficient as a function of x' where 
x is also fixed to 4. However in this case, a higher maxi- 
mum is found near the edge of the wave front on Fig. 
d). This suggests that the whole wave front is involved in the 
evolution. In fact, a better representation of the phenomena is 
found in the (x, x') plane of the sensitivity coefficient. On Fig. 
10, we present the contour plot of the sensitivity coefficient. 
The doted lines given on Fig. (0 correspond to the initial 
value of the density where the solution is perfectly correlated 
to the initial condition for Sh(x, 0)/6r)(x') — S(x — x'). For 
a given value of x' and coming from large values of x at time 
t = 2 (Fig. (13-a)), the sensitivity is still well correlated until 
one reaches a breaking off near the bottom of the wave front 
showing perturbations appear. After an interval, larger as the 
time elapses (t = 12, Fig. (J2j-b)), we observe a decorrelation 
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FIG. 2: Contour plot of the sensitivity coefficient 8h(x,i)/8r)(x') 
in the (x,x') plane, for two values of time (t — 2 and t = 12). 
For a smaller value of time (a), the solution is rather well correlated 
to the initial condition except in a small region near the wave front. 
When the time increases (b), the density spreads and one can observe 
a decorrelation through along the wave front. 



through along the wave front. In the mean time, the sensitiv- 
ity coefficient takes its larger values on this region. Behind the 
wave front, the density will be correlated again, although we 
constat a spreading when the time increases. For very large 
value of the time, the density is completely decorrelated as we 
said above with Eq. (II It . In order to sum up these results, one 
can say that the most the solution is sensitive to the initial con- 
dition, the most it is decorrelated in space, i.e. all the values of 
the initial condition participate to the solution at any distance 
along the wave front. This finally reveals a characteristic of 
the Burgers turbulence. 

Now we can find the sensitivity of the solution to the un- 
forced Burgers equation by the same method. In order to do 
the calculation, we have on one hand: 



8u(x, t) 
S(f(x') 



+ OC 



8u(x,t) 6rj(y) 
Sr)(y) Sip(x') 



dy = 



8u(x, f) 
8r)(y) 



(13) 



for we have Srj(y) / Sip(x') — 8(y — x' 
step function and on the other hand: 



i, where 9 stands for the 



8u(x, t) 
8r](y) 



8u(x, t) 8h(z, t) 
Sh(z,t) 8r)(y) 



dz 



-d x 



8h(x, t) 
Si](y) 



(14) 

These results allow us to calculate explicitly the sensitivity of 
the solution to the Burgers equation in term of Eq. JlOi . This 
is a rather cumbersome relation although simple to calculated, 
it is left to the reader. Moreover, the main conclusions about 
the sensitivity to initial conditions can be derived in the same 
manner as the one we have obtained above for h(x, t). 



IV. CONCLUSION 

The question arises now of the possibility to extend the re- 
sults issued from the preceding section to a larger class of 
equations. 

Let us consider first the forced Burgers equation by a pure 
time dependent term: 



d t u + ud x u = vd xx u + /(f). 



(15) 



Orlowski and Sobczyk fl2ll have found that an appropriate 
transformation onto the variable x and function u maps Eq. 
(II 51 to the unforced Burgers equation: 



y = x- (p(t), v(y, t) = u(x, t) - ip(t), 



(16) 



where 



<?(*)= / f (t) dr and ip(t) = [ <p(r)dr. (17) 

Thus the sensitivity coefficient of the h function in this trans- 
formation just undergoes a translation into Eq. dlOi : 



x ~ I f{i~)dT 
Jo 



(18) 



and does not affect the density otherwise. Then we are 
brought to the same conclusions as in the case of the unforced 
Burgers equation. 

Now we can give few additional remarks. The calculation 
of the sensitivity coefficient of Burgers and related equations 
may be quite naturally extended to the three dimensional case 
because the Hopf-Cole transformation still applies in this case, 
from which the densities may be easily deduced. There is also 
the case of initial/boundary problem for the Burgers equation. 
Actually, an analytic solution to this problem where solved 
several years ago and recently applied to related 

problems concerning the Burgers equation 1 16l ll7lfl8il . Sensi- 
tivity to initial conditions of such problems seems to be usefull 
to perform and will be the aim of a forthcoming paper. 
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